Failure of Scattering Interference in the Pseudogap State of Cuprate Superconductors 
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We calculate scattering interference patterns for various electronic states proposed for the pseudo- 
gap regime of the cuprate superconductors. The scattering interference models all produce patterns 
whose wavelength changes as a function of energy, in contradiction to the energy-independent wave- 
length seen by scanning tunneling microscopy (STM) experiments in the pseudogap state. This 
suggests that the patterns seen in STM local density of states measurements are not due to scatter- 
ing interference, but are rather the result of some form of ordering. 
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The origin of pseudogap phenomena is one of the 
key questions in the cuprate superconductors, as it pro- 
foundly affects their properties above the superconduct- 
ing transition, and potentially influences the cuprate 
phase diagram pi In a recent paper, we have used a scan- 
ning tunneling microscope to show that the electronic 
states in the pseudogap state of Bi^Sr-iC 'aC^Os+s are 
spatially modulated^ These electronic modulations in 
the local density of states, which are seen only for en- 
ergies less than the pseudogap energy scale, are ori- 
ented along the copper-oxide bond direction and have an 
energy-independent, incommensurate wavelength. One 
interpretation of these experiments is that these mod- 
ulations are caused by electronic ordering, a variety of 
which have been proposed for the cuprates in the pseu- 
dogap regime i^^Sii Another possibility is that such 
modulations are the consequence of scattering interfer- 
ence from one of the electronic states proposed for the 
pseudogap. 8 ' 9 The interest in scattering interference as a 
possible explanation is motivated by the success of this 
idea in understanding similar modulations in the super- 
conducting statei2*ii*i£, although even below T c there are 
a number of deviations from the scattering interference 
picturei 4 ' 1 ?! 14 ! 1 ^ In this paper, we elaborate on calcula- 



tions reported in Ref. 2 to demonstrate the failure of the 
scattering interference scenario in describing the modula- 
tions observed in the pseudogap state. We show that this 
failure is a generic feature of scattering interference itself, 
regardless of the model chosen for the pseudogap elec- 
tronic stateAiSiiLiSiiS This shortcoming suggests that 
the modulations observed in the pseudogap regime are 
the signature of some form of ordering. 

The scattering interference picture ascribes spatial 
modulations in the local density of states to interfer- 
ence between elastically scattered quasiparticlesiiS* 2 ^ A 
well-studied example of this phenomenon is the standing 
wave patterns seen on certain metal surfaces, in which 
electrons scattering from step edges and point impurities 
interfere coherently to form modulations in the density 
of statesi^i The period of these modulations at a given 
energy u> is determined by the wavevector q that joins 
the two points on curves of constant electron energy in k 
space with the greatest weight, as determined by the elec- 
tronic structure of the system. For a superconductor, the 
scattering interference contribution to the local density 
of states can be treated within the Born approximation, 
which introduces a correction, 
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5n(r,uj) = ~~^ Im { J d 2 fiG (r- f{,uj)V(fi,uj)G (fi - r, w) - F (f - f{, u))V(f[, ui)F (fl -r,u)} 
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to the density of states, where Go and Fq are the 
single particle and anamolous Green functions, re- 
spectively, and V is a weak, finite-range scattering 
potentiali 9 i 12 i 22 i 23 i 24 i 25 For pure potential scattering, the 
Fourier transform of the Born correction 6n(q, u>) — 
— -Im{V(q,ui)A(q, w)} seperates into a part 

A(q,u>) = J{d 2 kG {k,uj)G a {k + q,oj)- 

F (k,cj)F (k + q,ij)} (2) 
which contains all the wave interference information, and 
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d 2 xe- iq - x V{x,u) 
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which acts like a static structure factor^ For now, we 
will assume that the structure factor does not filter any 
wave interference information. 

In order to develop a context in which to understand 
quantum interference in the pseudogap state, we will first 
review the quantum interference picture for the two cases 
already discussed in the literature- the superconducting 
state and the Fermi liquid normal stateii 2 * 2 ^ 2 ^* 2 ^ Follow- 
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ing these previous works, scattering interference in the 
superconducting state can be modeled using the Green 
functions, 

^ n \ uj + iS + e,: 

G (k,uj) = - — - -it 

{u + i5-e % )(u + i5 + e % )-kl 

(oj + iS- e^)(uj + iS + eg) - Ai 

where eg = 120.5 — 595.1 x (cosfc x + cosfc a )/2 + 163.6 x 
cos k x cos k y — 51.9 x (cos 2k x + cos2fcj,)/2 — 111.7 x 
(cos 2k x cos fcy+cos k x cos 2fcj,)/2+51.0x (cos 2k x cos 2k y ) 
is the Bi2Sr 2 CaCu20s+s band structure from ARPES 26 
for slightly underdoped (S = .12) samples, Aj = 45.0 x 
(cos k x ~\~ cos/cy)/2 is the superconducting gap function, 
and 5 is a broadening term. As shown in Figures^ and 
b, this yields pictures with sharp peaks in q-space which 
disperse. To understand the origin of the dispersion qual- 
itatively, it helps to consider a simplified version of scat- 
tering interference, the so-called octet model, introduced 
by Hoffman et al. and McElroy et aZ.i 10 ! 11 Quasiparticles 
can elastically scatter between points on contours of con- 
stant electron energy in k-space (Fig. QJ;). The shape 
of these contours changes continuously as the energy 
increases. Consequently, the length of the characteris- 
tic scattering interference wavevector (q-space) changes, 
leading to dispersion. This picture changes dramatically 
if, instead of a superconducting Green function, we use 
a Fermi liquid normal state Green function by setting 
Ag = in Equation 4^ The sharp peaks in q-space 
for the superconducting state have been replaced with 
dispersing caustics (Figs. QJi and d). The equivalence 
of all points in k-space leads directly to an absence of 
sharp peaks in q-space. The dispersion is now a direct 
consequence of the band structure itself. Although the 
Fermi liquid scattering interference picture is useful in 
determining the effect of the band structure on scatter- 
ing interference, it is not applicable to the pseudogap, 
as it omits two key features of this state: the pseudo- 
gap in the density of states, and the ill-defined nature of 
quasiparticles. 

To extend the scattering interference model into the 
pseudogap regime, we must first choose an appropriate 
Green function with which to model the pseudogap state. 
Although its origin is still not understood, the pseudo- 
gap has been thoroughly characterized by ARPES in 
the underdoped cuprate superconductors £L At all ener- 
gies, peaks in the ARPES spectral function are so broad 
that it is unclear whether quasiparticles are still well- 
defined^. Meanwhile, at the Fermi energy, the spectral 
function shows extended arcs centered at the nodal points 
(Fig. |2K) in k-space. As energy increases towards the 
pseudogap energy scale, these arcs extend towards the 
Brilouin zone boundary. Theoretical attempts at under- 
standing the measured behaviour have either followed 
a phenomenological approach, or have proposed exotic 
electronic states for the pseudogap. We will first focus 
on the phenomenological attempts. Norman et al. have 




FIG. 1: (A) STM data in the superconducting state shows 
modulations in the local density of states which disperse. 
Shown here is data taken along the (0, n) direction (the Cu- 
O bond direction) from Vershinin et al~ at 40K for slightly 
underdoped Bi2Sr2CaCu20s+s (sqaures), and from Hoffman 
et aim- at 4K for underdoped (triangle), as-grown (diamond) 
and overdoped (circle) samples. The wavelength is in units 
of ao = 3.8A, the Cu-Cu distance. The solid black line shows 
the dispersion for the (0, tt) mode in the superconducting state 
calculated using Born scattering. The dashed black line shows 
the dispersion of the (0, n) mode for the Fermi liquid normal 
state. (B) The power spectrum of 5n(q, uf) in the supercon- 
ducting state with a broadening of lmV is shown here in the 
first Brillouin zone at an energy of uj = 20mV. Black regions 
correspond to a large contribution from scattering interfer- 
ence. The dispersion of the slowest dispersing feature along 
the (0, tt) direction, highlighted by the arrow, is plotted as 
the solid black line in (A). The Fq term in Eq. [I] makes 
the (tt,tt) mode have sharp peaks, but makes the (0,7r) mode 
weaker for non-magnetic impurities. The situation is reversed 
for the magnetic impurity case (not shown), in which the Fo 
term has the opposite signp^ (C) The solid black lines show 
contours of constant electron energy in the superconducting 
state at uj = OmV (point along (ir,n) direction), uj = A/2 
(small banana) and uj — A (large banana). The most promi- 
nent peaks in scattering wavevector space (q-space) along the 
(0, 7r) direction are indicated by the grey arrows. (D) The 
power spectrum of 5n(q, uS) in the Fermi liquid normal state 
with a broadening of lmV is shown here in the first Brillouin 
zone at an energy of uj = 20mV . The dispersion of the slow- 
est dispersing feature along the (0, n) direction, highlighted 
by the arrow, is plotted as the dashed black line in (A). 



modeled the ARPES data in the pseudogap regime using 
the Green function 

G (k,u) = (uj-ej:-^)- 1 (5) 
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and i^o = 0, where the self energy, 

S(fc,w) = -ir 1 + A|/(w + e^ + ir ), (6) 

Ti is a single particle scattering rate, T is a measure 
of decoherence, and A^ is a gap function with a k de- 
pendence which matches the Fermi arcsii& The calcu- 
lated scattering interference patterns (Fig. \2)p) have two 
notable shortcomings in comparison with the measured 
data shown in Figure^. First, the calculated patterns 
contain caustics, while the measured patterns contain 
discrete points along the (0, tt) directions. This short- 
coming can be overcome by assuming, for example, that 
the STM tunneling matrix element has a d-wave symme- 
try and selectively filters the density of states along the 
(0, 7r) direction. Second, the scattering interference pat- 
terns calculated from this phenomenological description 
still show features along the (0, tt) directions that dis- 
perse (Fig. [2Ji) - This dispersion curve resembles that in 
the superconducting state, as the characteristic wavevec- 
tor for scattering interference is of a similar length at 
the gap energy, but is shallower because this wavevector 
is much shorter at the Fermi energy in the pseudogap 
regime (Figs. ^ and|5^,). Even this shallower disper- 
sion results in a change of wavelength between 7ar> and 
4.6ao between the Fermi energy and the pseudogap en- 
ergy, whereas STM experiments show a fixed wavelength 
of 4.7 ± 0.2ao over this energy ranged 

We next consider the scattering interference scenario 
for various exotic electronic states proposed for the pseu- 
dogap regime. We will discuss how four of these exotic 
Green functions (two based on preformed pairs, and two 
based on ordering) produce dispersions which should be 
resolved in the data. Chen et al.— have proposed the 
Green function in Equation [3] with the self energy in 
Equation [5] to describe preformed pairs. However, unlike 
the Fermi arc model, they use a d-wave gap function. For 
small values of Tq, the scattering pictures resemble those 
calculated for the superconducting state, and for larger 
values, they resemble broader versions of those calculated 
for the Fermi liquid (not shown). Both these dispersions, 
similar to the ones shown in Figure la, should be resolved 
by the experiment. Another example of a preformed pairs 
calculation appears in Pereg-Barnea and FranaS, where 
the authors take a QED3 Green function Go(k,uj) — 
(oj + e^:) / (oj 2 + e| + Ai) 1 ^ 71 / 2 , with 77 being an anomolous 
dimension exponent which increases the amount of deco- 
herence of the patterns for larger values. As noted by the 
authors, the scattering interference patterns show a dis- 
persion identical to the superconducting state for i] < 0.5, 
which is shown here explicitly in Figures and b. We 
have also calculated the dispersion for a simplified d- 
density wave (DDW) ordering Green function Go(fc, oj) — 
(w + - -Afx— ejf )( w - A/x - eg +(3 ) - D 2 ), 

where D% = D (cosk x — cosfc a )/2 is the DDW gap, A/x 

is a chemical potential shift, and Q = (7r,7r)i2ii£ As can 
be seen in Figures and c, the simplified DDW Green 
function produces scattering interference patterns which 



disperse through a range of wavelengths (AA = 0.8ao be- 
tween -15mV and 35mV) which should be easily resolved 
by the STM experiment (maximum AA = 0.4ao over the 
same range of energies). Using a more complete DDW 
treatments^ results in an identical dispersion, and very 
similar scattering pictures (not shown). Finally, in Fig- 
ures |2k and d, we show the quantum intereference pat- 
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FIG. 2: (A) Shown in black are model representations of the 
electronic density of states in k-space measured by ARPES 
at the Fermi energy (arcs centered at the nodal points) and 
at the pseudogap energy (banana). Shown in gray are the 
shortest wavevectors along the (0, n) direction for scattering 
interference at the two respective energies. (B) We calcu- 
lated the power spectrum of Sn(q, w) for a phenomenological 
Fermi arc Green function where the extent of the Fermi arcs 
were chosen to give a scattering wavelength of 4.6ao at the 
Fermi energy, and unphysically small values of To = lmV 
and Fi = lmV were chosen to produce sharper features. The 
power spectrum of 5n(q, oj) for this Fermi arc picture is shown 
here in the first Brillouin zone at oj — 20m V . The dispersion 
of the most slowly dispersing feature along the (0, tt) direc- 
tion, indicated by the arrow, is shown in (D) as a solid line. 
(C) The power spectrum of the STM data at 15mV in the 
pseudogap state from Vershinin et al£ shows four peaks along 
the (0,7r) direction (circles). We have magnified the data to 
show only the relevant region in q-space. (D) These peaks 
in q-space correspond to a fixed wavelength of 4.7 ± 0.2ao, 
shown here as the grey line. Also shown is the dispersion of 
the most slowly dispersing feature along the (0, tt) direction 
in scattering interference calculations for the Fermi arc Green 
function. This feature disperses similar to the (0, tt) mode 
in the superconducting state for small (~ lmV) values of To 
and Fi. Larger (~ lOmV) values of Fo lead to the gradual 
suppression of this mode in favor of a mode which disperses 
similar to the band structure (Fig. ^). Larger values of Ti 
yield pictures which average over some part of the dispersion 
for a particular mode, leading to unphysically broad features. 
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FIG. 3: (A) The dispersions of all the exotic electronic state 
Green functions discussed in the text should be resolved by 
the STM experiment. Shown are the dispersions of the most 
slowly dispersing (0, 7r) modes for the QEDz (squares) , DDW 
(up triangles) and circulating current (down triangles) Green 
functions. These modes are labeled by the arrows in parts (B) , 
(C), and (D). STM data from the pseudogap state- is shown in 
gray. (B) The power spectrum of 5n(q, u) for a QED3 Green 
function is shown here in the first Brillouin zone at u — 20m V . 
The patterns were calculated using a value of rj — 0.4 and a 
broadening of lmV. For increased r\, the patterns in q-space 
become more diffuse until, by r\ ~ 1, there are no discernable 
features in q-space. (C) The power spectrum of 8n(q, ui) for 
a DDW Green function using a broadening of lmV and a 
chemical potential shift of +40meV is shown here in the first 
Brillouin zone at uj — 20mV. (D) The power spectrum of 
8n(q,Lj) for a circulating current Green function using A = 
.27, T = lOmV, D = WmV, and lu c = 500ml/ is shown here 
in the first Brillouin zone at w = 20ml/. Non-zero values 
of Do result in the slowest dispersing mode along the (0, n) 
direction (arrow) being replaced by a faster dispersing mode 
(circle) at energies uj < |2 x Do | . The dispersion shown in part 
(A) is for Do = 0, corresponding to the slowest dispersing 
mode. 



terns for the marginal Fermi liquid scenario with circulat- 
ing current order^S The Green Function for this phase, 



given by G (k,uj) = (ui-e n ± D(k) - Sj)" 1 for k ^ k f 
with £(fc,a->) = \{uj\ogx/Lu c + inx/ coshD(k)/x}, x — 
{uj 2 + tt 2 T 2 )(1/2) and D(k) = D (cos(k x a) - cos(k y a)) 2 , 
produces patterns which disperse in a fashion similar to 
the band structure. Thus, it also fails to match the dis- 
persionless feature seen in STM in the pseudogap regime 
(Fig. no. 

The failure of these approaches points to a fundamen- 
tal problem: any Green function will result in a wave 
contribution to Born scattering which disperses. The 
Green function is being asked to accomplish two con- 
tradictory tasks: on one hand, the imaginary part of the 
Green function (the density of states) must disperse in 
order to match the ARPES band structure, and on the 
other hand, the imaginary part of the Green function 
convolved with itself cannot disperse if it is to match 
the STM data in the pseudogap state. In order for Born 
scattering to simultaneously match the ARPES and STM 
data, the structure factor, assumed to be an all-pass fil- 
ter until now, could be chosen to pass only contributions 
near q = 27r/ao(0, 1/4.7) inside the pseudogap^ The real 
space scattering potential corresponding to such a struc- 
ture factor is an incommensurate, bond-oriented square 
lattice of scattering centers. While this can be justified 
on various physical grounds^&Si, it amounts to an ad hoc 
assumption that ordering exists. Moreover, any choice of 
Green function would then adequately describe the data, 
thus demonstrating that scattering interference fails to 
describe the interesting physics revealed in the experi- 
ment. 
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